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We study the relaxation processes of the infinitely long-range interaction model (the Husimi- 
Temperley model) near the spinodal point. We propose a unified finite-size scaling function near 
the spinodal point, including the metastable region, the spinodal point, and the unstable region. 
We explicitly adopt the Glauber dynamics, derive a master equation for the probability distribution 
of the total magnetization, and perform the so-called van Kampen Omega expansion (an expansion 
in terms of the inverse of the systems size), which leads to a Fokker-Planck equation. We analyze 
the scaling properties of the Fokker-Planck equation and confirm the obtained scaling plot by direct 
numerical solution of the original master equation, and by kinetic Monte Carlo simulation of the 
stochastic decay process. 



I. INTRODUCTION 



Relaxation phenomena in strongly interacting systems 
have been studied extensively, including various types of 
threshold phenomena. One of the most typical examples 
is the decomposition at the coercive field (the end of the 
hysteresis loop of a ferromagnetic system). This phe- 
nomenon appears in the field dependence of the order 
parameter in the ordered phase. In mean-field studies, 
this phenomenon is well described by the change in the 
free energy as a function of the magnetization (the order 
parameter of a ferromagnetic system). That is, in the or- 
dered state, the free energy has two minima representing 
the symmetry-broken states. When we apply an external 
field, one of them is selected to be the equilibrium state. 
For a weak field, however, the other state remains as a 
metastable state. When the field becomes strong, the 
metastable state finally becomes unstable. This point 
is called the spinodal point. However, in models with 
short-range interactions, fluctuations cause the system to 
escape from the metastable state through localized nucle- 
ation phenomena, and thus the change of the relaxation 
time is smeared. Thus, although there are several stud- 
ies on the divergence of the relaxation time near spinodal 
points 111, the spinodal phenomenon in short-range sys- 
tems must be defined only as a crossover [2j]. 

Recently, however, it has been pointed out that the 
mean-field universality class is realized in spin-crossover 
type systems ^3|. In these materials, the volume of the 
unit cell changes, depending on the local two-fold states 
(say the high-spin state and low-spin state). The vol- 
ume change causes a lattice distortion, and the elastic 
interactions among the local distortions cause an effective 
long-range interaction among the spin states. The criti- 
cal properties of the spin-crossover system were found to 
belong to the mean-field universality class. It was also 
found that the finite-size analyses for various quantities 
are very similar to those of the long-range, weakly inter- 



acting model (the Husimi- Temper ley model). 



In spin-crossover systems, one may expect that the dy- 
namics, as well as the equilibrium properties, corresponds 
to that of the mean-field model. In particular, various 
threshold phenomena have been pointed out in the dy- 
namics of the spin-crossover type materials. For example, 
a spinodal phenomenon without nucleation type cluster- 
ing was reported in numerical calculations [3|, The 
change is not a crossover process, but rather a change 
with a true critical singularity that can be described by 
mean- field dynamics 5]. The metastable state does not 
relax to the stable state in the mean-field approxima- 
tion, which corresponds to the infinite system size in the 
Husimi- Temperley model. In recent extensive studies on 
nano-size systems, finite-size effects turn out to play im- 
portant roles. Therefore, it would be interesting to study 
finite-size effects on the spinodal phenomenon as a critical 
dynamical process. The phase transition near the mean- 
field spinodal is well known [Ij, and a numerical study 
has been reported recently 0. However, a unified scaling 
function for the spinodal lifetime has not yet been explic- 
itly given as far as we know. Therefore, we here study it 
in the long-range model with standard Glauber dynam- 
ics. We first derive a master equation as a function of 
the total magnetization, which is possible in the Husimi- 
Temperley model because of the long-range nature of the 
interactions. Then, we derive a Fokker-Planck equation 
by using an expansion in terms of the inverse system 
size, which is an example of the van Kampen Omega- 
expansion (7i, iSj] . We analyze the Fokker-Planck equation 
and derive a scaling relation and also asymptotic forms 
of the relaxation times. These properties are confirmed 
by direct numerical investigations of the original master 
equation, as well as corresponding Monte Carlo simula- 
tions of the stochastic decay process. 
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II. INFINITE LONG-RANGE MODEL AND 
THE SPINODAL POINT 

We investigate the relaxation phenomena near the 
spinodal point in the infinitely long-range model. This 
is a spin model in which each spin interacts equally with 
every other spin. The Hamiltonian is 



Ti = — —M'^ - HM, M = V 



N 



2N 



(1) 



where H is the magnetic field and (Ti ~ ±1. It is well 
known that the mean-field theory is exact for this model 
in the limit of iV — > oo. This model shows a second- 
order phase transition at T = Tq = J, H = 0. Below the 
critical temperature, a metastable state exists for weak 
magnetic fields. When the magnetic field becomes strong, 
the metastable state becomes unstable at a certain point, 
known as the spinodal point. In order to determine the 
spinodal point, we consider the extended free energy, i.e., 
the free energy for fixed total magnetization. In the in- 
finitely long-range model, this is given by 

f{m) = -^"^^ - - ^lnwC(Ar+M)/2 

^ — Hm 

2 

1/l + m 1 + m 1 — TO 1 — m\ 

+ - n —— + n , (2) 

f3\ 2 2 2 2 J ' ^ ' 

in the limit N oo where we use Stirling's formula 
and m is the magnetization per spin (m = M/N). The 
spinodal point is given by the following conditions. 



_L = o and ^=0, 



which give 



1 



Hsp = - — ± — In 

2/3 i_ /i_ 1 



f3J 



at which the magnetization is given by 



msp 



±A 1 



1 



(3) 



(4) 



In this paper, we consider the case where we increase the 
field from a negative value. Therefore, we consider the 
behavior of a locally stable state at negative magnetiza- 
tion around tosp < 0. 



III. THE DYNAMICS NEAR THE SPINODAL 
POINT 

We study the dynamics via the standard master equa- 
tion 

^E^^-J2ws^s'PiS) + Y,Ws'^sP{S'), (5) 



where 5* and S" denote states of the system and Ws^s' 
is a transition probability from S to S' . The probability 
of the state S at time t is denoted by P{S, t). 

Among the many possible dynamical models (choices 
of the transition probability), we adopt the Glauber dy- 
namics in this work. In the Glauber dynamics, the transi- 
tion takes place as a flip of a local spin, and the transition 
rate Wij from a local spin state i to a local spin state j 
is given by 



Wi. 



To e 



-f3E, ' 



(6) 



where Ei denotes the energy of the system in spin state 
i, and tq is some characteristic time scale. In this paper, 
we scale the time by tq and set tq = 1 for simplicity. 

With this transition rate, we construct a master equa- 
tion for the mean-field model. As the Hamiltonian de- 
pends only on the magnetization M, the master equation 
is written in closed form for M. Thus the master equa- 
tion ([5]) can be expressed as a function of M. Let P{M) 
be the probability that the system has the total mag- 
netization M . The master equation for P{M) is given 
by 



9P(M, t) 



1 



dt To 

N + M exp[-/3( J(M 



l)/N + H)] 



2 2cosh[l3{J{M -l)/N + H)] 
N-M (iyLY>[l3{J{M + l)/N + H)] 



P{M) 



2 2cosh[l3{J{M + l)/N + H)] 
N~M + 2 exp[/3(J(A/- l)/iV + i/)] 



P{M) 



N- 



2 2cosh[(3{J{M ~1)/N + H)] 
M + 2 exp[-/3(J(Af + 1)/A^ + H)] 



P{M 



2co8h[P{J{M + l)/N + H)] 



P{M 



'2) 
2) 



(7) 



where N is the number of spins and M takes discrete 
values -TV, -N + 2,---N (see Appendix E|. For fi- 
nite N we can solve the simultaneous equations for 
P(-iV, t), P{-N + 2,t),--- , P{N, t), as well as perform 
a Monte Carlo simulation of the model. 



IV. ANALYSIS OF RELAXATION NEAR THE 
SPINODAL POINT 

Next, we study the scaling properties of the relaxation 
time near the spinodal point. For large N, we put to = 
M/N and regard m as a continuous variable. Expanding 
the RHS of Eq. ([7]) in a series in e = we obtain the 
following Fokker-Planck equation: 



dP{m) 
dt 



d d'^ 
— 5i(to)P(to) -f- £— ^52(m)P(m) + O(e^) 
am om'' 



(8) 
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where 

gi{m) = m — tanh[/3( Jm + H)] + s 
32 (™) = 1 — mtanh [f3{Jm + H)] . 



(3Jm 



cosh^[/?(Jm + i/)] 



(9) 



The last term of gi (m) gives the correction to the spin- 
odal point due to the finite-size effect. Hereafter, we ig- 
nore this term as it is very smalL Near the spinodal point, 
we expand gi(m) and (72 (w) around the spinodal point. 
We set X = m — msp, and y — (3{H — i?sp) ^ h — hgp, 
and then we have 



gi{m) ~ ~jj ~ rjipjx^ + 2xy) 
52M ^ i^. 



(10) 



where = |tosp| = yjl — \/ (iJ. 

Let us consider the time evohition of x, starting from 
X = Q. The distribution of x evolves according to Eq. ([8|) 
with gi ~ —y/fiJ and 32 ~ 1/(3 J. When x approaches 
X ^ y^^'^, the correction term —r]PJx^ in gi becomes 
relevant, but other correction terms in gi, such as —2rjxy, 
are still irrelevant. Therefore, in order to determine the 
relaxation time, we can use the approximation that gi ~ 
—y/ j3J — rjPJx^ in the early stage of the phase change. 
Then, the Fokker-Planck equation takes the form 

(11) 

Now, we introduce the scaled parameters, 

^^Ay\-"\ (12) 



and 



A = ^2/3 



y- 



(13) 



Then, for nonzero y, the Fokker-Planck equation is given 

by 



1 



(14) 



where A > for the upper sign, and A < for the lower 
sign. Because Eq. (fT4|) depends only on A except for the 
factor e^^^, the relaxation time is expected to be given 
in the form 

T ^ 7Vi/V(A) = iVi/3/(iV2/3(/i - hsp)). (15) 

This is the finite-size scaling for the relaxation time near 
the spinodal point. 

However it is necessary to pay attention to the range 
of the parameters, in which the application of the above 



estimate can be verified. We may regard the relaxation 
time as a time when the magnetization x becomes 0(1). 
Then, ^ becomes 0(|y|~^/^). This implies that we may 
need to include an additional y-dependence in the relax- 
ation time. Thus, wc cannot immediately conclude that 
the finite-size scaling has the simple form of Eq. p5|) . 
This problem is investigated in Appendix |b1 where we 
show that the additional contribution does not need to 
be taken into account, and the finite-size scaling is cor- 
rectly given by Eq. (|15p . In the following, we investigate 
the relaxation time for the cases y < 0, y = 0, and y > 0- 



y < case 




FIG. 1: The scaled potential U{^) in Eq. dTT)). A: the 
metastable state. B: a point in the basin of attraction of 
the stable state. C: the unstable state. 



This case corresponds to the relaxation from the 
metastable state to the equilibrium state. We can esti- 
mate the relaxation time by Kramers' method of escape 
over a potential barrier 0. We rewrite Eq. in the 
following form (A < 0) 



dP^ 1/3^ 

dt ^ di \(3J\K\ 



1 /3J|A|3/ = ;7(C)^g/3J|A|-''/^C/(Op 



where 



(16) 



(17) 



is the scaled potential. We consider the escape of proba- 
bility from the valley to the outside (A — > B) as depicted 
in Fig.[T] 9]. There, the probability current a is given by 



We consider the stationary current [a =const.) and in- 
tegrate it between the two points A and B. 



£l/3 [e/3J|A|^/=C/(0p]-B 



(19) 



If the number of spins N is very large, we can use the 
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steepest-descent method, and we have 

i-B 



„/3J|Al^/^C/(C) 



„/3J|A|3/2(7(C) + i/3,/|A|=*/=;7"(C)(C-C)=' 



'2-K 



/3J|A|3/2f7"(C) 



„/3J|A|^''^;7(C) 



(20) 



This approximation is vahd for a sufficiently large |A|. 
We consider the early stage of the relaxation and assume 
the relaxation has not occurred yet, namely P{B) sa 0. 
Then, we obtain the estimate: 



a = P{A)e^ 



/3 / U"{C) pj\\\^/^(u(c)-U(A)) 

27r/3J|A|i/2 



(21) 



The probability distribution near the point A is approx- 
imately given by 



P{0 ^ nA 



exp [-/3J|A|3/2 |[;(^) ^ iC/"(A)(C - A)2}] 



where Z is a partition function, 

/>OC 

Z = I exp 



(22) 



;3J|Ap|c/(A) + ic/"(A)(C-A)^ 

(23) 

Namely, the probability distribution near the point A is 
given by the equilibrium distribution of the approximate 
harmonic potential. The variable ua represents the total 
probability near the point A. This quantity evolves as 



d 

dt""^ 



1 



(24) 



where r is the relaxation time that we want to know. As 
the probability at the point A is given by 



P{A) ~ riA 



exp [-/3J|A|3/2[/(A)] 



Z 



HA 



pj\A\y^U"{A) 
2tt 



(25) 



(see Eq.Jlll)), combining Eqs.(l2T|), dH]), and ([25]), we ob- 
tain 



- N^/^2tt\U"{A)U"{C)\-^/^\A\-^ 
xexp [/3J|A|3/2(t/(C) -C/(A)) 



-1/2 



(26) 



This is the result of the well-known Kramers' formula for 
the escape rate, and it agrees with the finite-size scaling 
equation pS]) . 

The potential U{0 is C/(0 = UPJ-vPJ^V^, and the 
two points A and C are given by the condition dU /d^ = 0. 



Therefore, 



A 



-C 



/5 Jryi/a 



U{C) - U{A) ^ \ ^ 



3 (/3J)27,l/2 ' 

V"{C) - -2??i/2. 



(27) 



Hence, the relaxation time for sufficiently large |A| is 



(28) 



Here, it should be noted that from Eq. ^ the following 
relation holds: 

A3/2 = = p{F{C) - F{A)). 



Therefore, the relaxation time obtained above is roughly 
T ^ e^*^^, which is the well-known Arrhenius formula. 
Another derivation of Eq. ([28| uses the WKB approxi- 
mation as discussed by Tomita, et al. ■ We can derive 
the same result (see Appendix [Ct . 

In the infinite long-range model, microscopic fluctua- 
tions do not grow to become macroscopic because the 
long-range interaction suppresses clustering and nucle- 
ation. It should be noted that, in the limit of iV ^ oo, 
the system remains at the metastable or marginally sta- 
ble point, and the relaxation time from that point be- 
comes infinite. In contrast, in systems with short-range 
interactions, the nucleation process causes the system to 
relax to the equilibrium state in a finite relaxation time. 
Thus, the divergence of the relaxation time does not take 
place, and so far the divergence of the relaxation time 
has not been considered seriously. However, as has been 
pointed out effective long-range interactions appear 
in systems in which elastic deformation mediates interac- 
tions among the spins. In such systems, the long-range 
interaction model is effectively realized, and we expect 
that the finite-size scaling discussed here would be rele- 
vant. 



y = case 

Next, we consider the relaxation just at the spinodal 
point, y = 0. Substituting y = in the Fokker-Planck 
equation (fTTj) . we obtain 

|P(.,0=./34(.^P(.,^))+|^^P(.,0. (29) 
Putting X ~ e^^'^z, 
|p(.,0 = .^/-^{-./3j|.^ + ^^}p(.,0. (30) 

By using the scaled variable s = te^^^, we can eliminate 
the e-dependence. Thus, as pointed out by Kubo et al. 



5 



1 , the relaxation time behaves as 



(31) 



The relaxation time diverges in the limit of iV ^ oo 
just at the spinodal point. In the limit of — > cxd, the 
system remains at the unstable point, and the relaxation 
time becomes infinite. This divergence is again due to 
the long-range interaction. 



2/ > case 

Finally, we consider the case y > 0. In this case, even 
\i N = oo (e = 0), the relaxation takes place. Namely, 
the relaxation time saturates at a finite value at large 
N . Therefore, we consider only the limit N oo. The 
Fokker-Planck equation (|14p then becomes 



dP 
'dt 



-y 



1/2 



d 



1 



(32) 



dS, \(3J 

Therefore, the relaxation time is expected to scale as 

r^y-i/2~(/i-/jSp)-i/2. (33) 

In the limit of — > cx), there is no diffusion term, so we 
can derive the time evolution of the scaled magnetization 
C(t) directly. Namely, putting P(C,t) = S{C - ^{t)), we 
obtain 



'im(t)'- 



The solution of Eq. (l34l) is given by 

1 



tan 



riyt 



(34) 



(35) 



for yjrjyt < 7r/2. At a time \/r]yt — 7r/2, the above ex- 
pression indicates that ^(t) would diverge. However, the 
higher order terms in the original Fokker-Planck equation 
([5]) prevent this divergence of the magnetization. Thus, 
the relaxation time is estimated as 



In the scaling form, 



7Vi/3|(^A)-i/2 



(36) 



(37) 



This result is consistent with that obtained by Binder [l[ . 
He showed that for h > /isP: the relaxation time behaves 
asT~ (/i-/isp)"^/^- 



V. NUMERICAL RESULTS 

We have seen that the relaxation time r obeys the 
finite-size scaling form (|15p : 



r{v,N)=N^'^f{N^l^v), 



and we have derived the asymptotic forms of the relax- 
ation time, i.e., 



/(A) 



^(?7|A|) 



-1/2 



for A > 1 



3/2 



for - A > 1 



(38) 

We checked these results by solving the original master 
equation ([7]) numerically, and also by performing kinetic 
MC simulations. The parameters are set as /3 = 1 and 
J = 2. The relaxation time is defined as the time at 
which the magnetization of a sample reaches a certain 
value mg. Here we adopt toq = 0- In the Monte Carlo 
simulations, the relaxation time is measured directly in 
each realization. On the other hand, in the master equa- 
tion we have to define it from the change of the probabil- 
ity distribution P{M,t). Namely, we obtain the average 
of the relaxation time with the formula: 



dt 



M<0 



P{M, t)t. 



(39) 



We plot data for various parameters in a scaling plot 
in Fig. [2l namely, in the coordinates (A = N^^^y, 
ln(A^-i/'^r)). We obtained data by solving the master 
eqaution using the above formula, and also by perform- 
ing corresponding kinetic Monte Carlo simulations (see 
Appendix D). We confirmed that both methods give the 
same results, as they should. In the Monte Carlo simu- 
lations, each data point is an average over 1000 samples, 
and the error bars are smaller than the symbol size in the 
figure. All the data collapse well onto a scaling function, 
which indicates that the finite-size scaling works well. 
For large A, data points for different N deviate from the 
scaling function. This fact is explained as follows: The 
condition for the finite-size scaling to hold is that the 
system size N is sufficiently large and the magnetic field 
is sufficiently close to the spinodal point. This implies 
A^ ^ 1 and \y\ <g; 1. However, an even stronger condi- 
tion is required for the finite-size scaling. As we assumed 
X = m — msp <C 1 to derive the Fokker Planck equation 
and X was rescaled as x = CIj/T^^i not only \y\ <C 1 but 
also |y|i/^ ^ 1 was necessary. Therefore, 

|A|l/2 ^^1/3 

is required for the finite-size scaling. 

Now, we compare the numerical results of the master 
equation with the asymptotic form for the relaxation time 
r, Eq. l(38|) . In Fig. [2l we compare numerical results 
and the asymptotic forms, Eq. ([38| . Here we find that 
the asymptotic formulae describe the scaling form well in 
the large-A region, where the data points appraoch the 
asymptotic forms when the size increases. Here, we find 
that the scaling property (fT5|) holds, and the asymptotic 
forms also hold asymptotically. 
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FIG. 2: Comparison of the asymptotic form of the relaxation 
time with numerical results. Data are plotted in the coordi- 
nates (A, ln(7V"^''^r)). ME and MC represent the solutions 
of the master equation and the Monte Carlo results, respec- 
tively. The solid lines denote the asymptotic form 



VI. SUMMARY AND DISCUSSION 

Mean-field type critical behavior takes place in sys- 
tems with effective long-range interactions, as has been 
pointed out for spin-crossover type materials [sl] . We ex- 
pect that the dynamical critical properties such as the 
spinodal phenomena are realized in those systems. In 
models with short-range interactions, there exists a mode 
of relaxation from the metastable state through nucle- 
ation of localized clusters. Thus, the relaxation time 
around the spinodal point changes smoothly, and the 
critical properties at the spinodal point in the mean-field 
theory are smeared out. However, in long-range inter- 
action models, the relaxation time diverges as describled 
by the mean-field theory. It is, therefore, necessary to 
study the finite-size scaling properties of the critical be- 
havior. Thus, we here studied the size dependence of the 
relaxation time near the spinodal point in the Husimi- 
Temperley model. We derived the master equation for 
the probability density of the total magnetization under 
the Glauber dynamics, and from it we derived the Fokker- 
Planck equation by the Omega-expansion method. Using 
this Fokker-Planck equation, we investigated the relax- 
ation processes near the spinodal point. As a result, we 
obtained a finite-size scaling function for the relaxation 
time, which covers both sides of the spinodal point, i.e., 
the metastable side and the unstable side. We also de- 
termined the asymptotic forms of the scaling function. 

The critical properties obtained in the preset work 
should apply widely to threshold phenomena in long- 
range interacting models, such as the threshold phenom- 
ena found in the excitation process by photo-irradiation 
from the low temperature phase to a photo-excited high- 
temperature phase in spin-crossover materials We 
hope the scaling properties presented here will help to 
analyze such processes in experimental systems. 



APPENDIX A: THE DERIVATION OF THE 
MASTER EQUATION 

In this Appendix, we derive the master equation ([71) for 
the Hamiltonian ([T]) and the transition probability 
The probability of the state {cti, (72, • • • , cttv} at a time i, 
which is denoted by P((Ti, • ■ • ,crAr;i), evolves according 
to 



d_ 

dt 



N 



P{ai, ■ ■ ■ ,aN;t) = -^WA/((Ti -(Ti)P(o-i, • • • ,(TN;t) 



N 



+ ^ C^Af _2(Ti {-CTi <yi)P{- ■ ■ , -O-i, • • • ]t) (AI) 



where 



1 exp (-/3(T,(J(M - a,)/N + H)) 



TO 2cosh{(3{J{M - a,)/N + H)) ' 

(A2) 

and M is the total magnetization, i.e., M — J^i'^i- We 
consider the time evolution of the probability of M : 



P{M, t) 



J2 <5(5]a„AnP((Ti, 



1 (^N;t), 



where 6{a, b) denotes the Kronecker delta. After some 
calculation from Eq. (jAip . the equation of motion for 
P{M,t) is obtained in the form 

d N + M 



N - M 



LOMi-l^+l)P{M,t) 



N -{M - 2) 
2 

N + (M + 2) 



^M-2(-l^+l)P(Af-2,t) 

ujM+2{+l ^ -l)P{M + 2,t). 

(A3) 



The meaning of this equation is clear. The first and sec- 
ond terms correspond to the transition from the magne- 
tization M to M — 2 and M + 2 respectively. The third 
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and fourth terms represent the transitions from M — 2 
to M and from M + 2 to M, respectively. This equation 
gives Eq. dl]). 



APPENDIX B: CUT-OFF INDEPENDENCE OF 
THE FOKKER-PLANCK EQUATION (fT4l) 



In Sec. Illlj we remarked on the possibihty of an ad- 
ditional y-dependence in the relaxation time. Namely, 
if we regard the relaxation time as the time when the 
magnetization x becomes 0(1), this corresponds to the 
time when ^ becomes 0(|i/|~^/^), and this implies that we 
cannot conclude the finite-size scaling of the relaxation 
time, Eq. ([15]), from the form of the Fokker-Planck equa- 
tion p^ . In other words, although the Fokker-Planck 
equation ([T4|) seems to depend only on A, we must re- 
strict the range of the variable ^ < 0{y^^^^) and this 
cut-off oi ^ can induce an additional y-dependence in the 
relaxation time. We note that the relaxation time indeed 
depends on the cut-off in other situations. One exam- 
ple is the relaxation from the mean-field unstable fixed 
point. In this case, the Fokker-Planck equation is given 
by 



d_ 

dt 



P{x,t) = -^xP{x,t)+e^P{x,t), (Bl) 



(we set some coefficients equal to unity). We can trans- 
form this equation to the scaling form similarly. If we set 



d_ 

dt 



d_ 
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9e 



(B2) 



This equation is apparently independent of e. Is the re- 
laxation time independent of the system size N — 1/e? 
The answer is No. It is known that the relaxation time 
in this case is t ~ In TV [ll|. We show that this TV- 



dependence stems from the finite cut-off. We can solve 
Eq. im for the initial condition P(C,0) = (5(C), 



where a{t) is given by 



1 



27rCT(t) 



cxp 



2a{ty 



a{t) = ^fe 



2t 



(B3) 



(B4) 



It takes infinite time for ^ to reach infinity. As x = e^^^C, 
x'^{t) ~ £a{tY = ee^*. We consider the relaxation time 
as the time when x'^{t) reaches 1, i.e. x'^{t) ^ 1, the 
relaxation time is proportional to In TV, 

r-lnTV. (B5) 

In this way, we found out that the cut-off dependence 
could actually affect the relaxation time, but this cut- 
off played no role in the case of the relaxation near the 
spinodal point. 



Hence, here we show that the relaxation time does not 
depend on the cut-off of ^ if this cut-off is very large. 

If we denote the average of C" over P{S,)hy C", the time 
evolution of ^ is given by 



^(t)=ei/3|A|i/2(± 



> e 



l/3|A|l/2 



± 



1 

1 

Jj 



vP-mf]- (B6) 



If £,{to) is larger than 1/77^/^/3 J, it can be shown from 
Eq. jTO]) that 



1 - 



|ggexp(2asV3AV2(t-^o)) ^ 



exp(2a£i/3Ai/2(t -to)) 



(B7) 



where a — The average of the scaled magneti- 

zation C(t) reaches infinity when the denominator of the 
RHS of the above equation is zero, namely 



t-tn 



2,i/3|A|i/2 V?(io)-a 



(B8) 



Because <o is finite, it takes only a finite time for ^{t) 
to reach infinity. Therefore, there is no cut-off depen- 
dence on the relaxation time in the Fokker-Planck equa- 
tion (HH). 



APPENDIX C: DERIVATION OF EQ. ([MD BY 
THE WKB APPROXIMATION 




A 

metastable 



equilibrium 




equilibrium 



FIG. 3: Rough sketches of the scaled potential U{^) and the 
corresponding Schrodinger potential V{^). 



In the body of this paper, we estimated the relaxation 
time for A < and |A| ^ 1 according to Kramers' ar- 
gument. Here we give another derivation by using the 
WKB approximation. The following derivation is es- 
sentially the same as that of Tomita, et al. The 
Fokker-Planck equation ^B)) can be transformed to the 
"Schrodinger equation" 



rl/3 Q2 



dQ ^ 

dt ~ pj\A\ de 



^^'^V{OQ{0 ^ e^'^HQiO (CI) 
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by substituting 

P = exp (^-i/3J|A|3/2[/(o) QiO- 

The scaled potential C/(^) is given by Eq. (|17p . and the 
Schrodinger potential V{£_) is 

ViO = \pj\A\^U'{0' - \\A\'/^U"{0. (C2) 

If the eigenvalues of TC are (i = 0, 1, 2, • • ■ ) and the 
eigenfunctions are 0i, we can expand Q{S.,t) as 



Q(C,i)-^c,(/.,(e)e- 



(C3) 



The lowest eigenvalue is Aq = 0, and the corresponding 
eigenfunction is 



MO 



1 



Zl/2 



(C4) 



which corresponds to the equilibrium state. Z is the 
normalization factor for J (pQd^ = 1. The second lowest 
eigenfunction will represent the metastable mode, and 
the corresponding eigenvalue will be connected with the 
inverse of the lifetime of the metastable state, e^/^Ai ~ 

From the Schrodinger equation (jCip . 



= -^<J>'^ + V<l)o 
2m 



Ai0i 



2m 



01 +V^0i. 



(C5) 
(C6) 



The mass m is m = /3J|A|/2. If we multiply Eq. (jC5[) by 
01 and Eq. (|C6[) by 0o and subtract the two equations, 
we obtain 



Ai0o0i = 2m ^'^"'^^ ~ 01 0o)' 



(C7) 



Integrating this equation from — cxd to the point C (see 
Fig. [3]), we obtain 



Ai 



(C8) 



because 4>'q{C) = 0. Here we consider the metastable 
wave function 0nis, which corresponds to the localized 
canonical distribution at the valley of the potential. 
Hence we assume for ^ < C 



00 ~ w0ms for ^ < C, 
where m is a constant given by 



c 



pj\A\y^U"{A) 



(C9) 



(CIO) 

Besides we assume that the first excited eigenfunction is 
also proportional to 0ms in the range £. C, 



t>i « vq 



for e<C, 



(CU) 



because 0i is considered to represent the metastable 
mode. Under these assumptions, we get 



0O0ldC = UV 0ms ~ UV 



(C12) 



Using the WKB approximation, it is obtained that 



0'i(C) « ^2mViC)MC) = -V2mU(C)0o(C). (C13) 



Substituting Eqs. ((UT2)) and (|UT3)) into Eq. (|U8|) . 



V 2m 

After some calculation, we obtain 



(C14) 



Ai - y^|C/"(A)[/"(C)|e-'5^^ (C15) 

where the free energy barrier /3AF is fi/^F = 
/3J|A|3/2(f7(c) _ U{A)). Therefore the lifetime of the 
metastable state r, which is equivalent with the relax- 
ation time, is 



T ~ iVi/3y2^|C/"(A)t/"(C)|-i/2|A|-i/2e/3AF_ (^16) 

Comparing with the relaxation time obtained by 
Kramers' argument, Eq. ((26|) . they agree with each other 
except for the minor difference in the constant prefactor. 



APPENDIX D: MONTE CARLO SIMULATION 

We also obtained data by performing kinetic Monte 
Carlo simulations to confirm the data obtained by solving 
the master equation. In the Monte Carlo simulations, 
each data point is an average over 1000 samples, and the 
error bars are smaller than the symbol size in Fig. [21 

The algorithm of the Monte Carlo simulations is as 
follows. We choose a spin at a site i randomly, and up- 
date the spin with the probability corresponding to the 
Glauber model given by Eq.®: 



exp {-l3a,{J{M - a,)/N + H)) 
2 cosh J(M ~a,)/N + H)) 



X At, 



where <Ti is the spin state of the i-th spin (ai ~ ±1)- In 
principle, a small time increment At <C 1 is neccesary 
to reproduce the result of the master equation given as a 
differential equation. However, we found almost the same 
result with the different time division, At=0.01 and 1 for 
the quantities plotted in Fig. [21 Hence we obtained the 
data with At = 1 . During a Monte Carlo step we perform 
single spin flips N times. Therefore the time t is related 
to the Monte Carlo step s by t = s because At = 1. The 
initial condition of each Monte Carlo simulation is set to 
the spinodal magnetization. 
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